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Abstract We present a multi-frequency and multi-instrument study of the 20 
January 2005 event. We focus mainly on the complex radio signatures and their 
association with the active phenomena taking place: flares, CMEs, particle accel- 
eration and magnetic restructuring. As a variety of energetic particle accelerators 
and sources of radio bursts are present, in the flare-ejecta combination, we 
investigate their relative importance in the progress of this event. The dynamic 
spectra of Artemis-IV-Wind/ Waves-Hiras with 2000 MHz-20 kHz frequency 
coverage, were used to track the evolution of the event from the low corona to 
the interplanetary space; these were supplemented with SXR, HXR and 7-ray 
recordings. The observations were compared with the expected radio signatures 
and energetic-particle populations envisaged by the Standard Flare-CME model 
and the reconnection outflow termination shock model. A proper combination of 
these mechanisms seems to provide an adequate model for the interpretation of 
the observational data. 
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1. Introduction 

Solar flares and coronal mass ejections (CMEs) are the most energetic phenom- 
ena on the Sun. They can be envisaged as t wo different aspec ts of a common 
magnetic energy instability and release (e.g. IPick et al. . 2006 . and references 
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therein), described by the standard CME- flare model which combines the kine- 
matics of the energy release of ejecta and flares: As the rising CME (or erupting 
filament or rising plasmoid) stretches and deforms the coronal magnetic-field, 
a vertical current sheet is formed, where recon nection explosive l y releases the 
free magnetic energy, stored i n the corona (e.g. Forbes and Lin . 2000t Forbesl . 
2003j iPriest and Forbesl . I2002T) ; furthermore, CME accele ration profile and flar e 



ene7gy release are seen to^olve in a synchronized manner femmer et oilM 



The energy, thus liberated, is divided over plasma heating, particle acceleration, 
kinetic energy of the eruption, and MHD shock waves. 

The development, in the low corona, of large flare/CME events coincides 
with an extended opening of the magnetic field, accompanied by energetic- 
particle acceleration and injection into interplanetary space and shocks. These 
ar e detectable, at metric a nd longer waves, by their radio signatures (e.g. review 
bv lPick and Vilmerl . l2008h . 

The Type III burst radio emission is produced by non-thermal electrons 
streaming along coronal magnetic lines. In open field li nes, they may escape 
into interplanetary space and may be detected in situ (e.g. | Faindberg and Stond. 



1970 : Lin. Evans, and Fainbere . 1973 ; Kurt. Logachev. and Pissarenkol 19771 : Klein et aZj j 



20081 ). In closed magnetic structures, on the other hand, they eventually turn 
Sunwards resulting in inverted U- or J-shaped bursts on the dynamic spectra (U 
or J bursts) . Metric bu rsts of the Type III family generally consist of groups of 



individual bursts (e.g. iGoldmanl 11983. and references within); on the reports 



they are marked as III G for fewer than 10 occurences and III GG for more. 

The Type II bursts trace the propagation of MHD shocks in the corona 
and interplanetary space. It is, in general, accepted that Type II bursts at 
decametre and longer w avelengths are driven by CMEs, bow or flank, (e.g. 
Vrsnak and CliverL I2008L for a review). The Type II bursts at metre wavelengths, 
howe ver, are interpreted eithe r as a flare blast wav e ( Vrsnak. Magdalenic. and Aura ss.l 
200lh or a CME driven shock (|Kahler et adll984lMaia et aUbOOOUClafien and Aurassi | 

2002). A metric stationary Type II burst (a Type II with a very slow drift on the 

dynam ic spectrum) has been identified, for the first time, bv lAurass. Vrsnak. and Mann^ 
(2002); it was interpre t ed as the radio signature of a fast- mode termination 
shocks bv lAurass et al. ( 2006| ). The proposed formation process included a pair 
of counterstreaming outflow jets moving upwards and downwards from the recon- 
nection site, encountering the rear of the C ME and the top of the post flare loop, 
respectively, resulting in termination shocks Jm ann. Aurass. and Warmuthl . [2006:1 
Mann. Warmuth. and Aurasst l2009t IWarmuth. Mann, and Aurassl l2009h . 

The continua during periods of activity, represent the radiation of energetic 
electrons trapped within magnetic structures and plasmoi ds, and they are know n 
under the names of Type IV bursts and "Flare Continua" (|RobinsonLll978[ll985h . 
The stationary Type IV (IV mB) bursts emanate from magnetic structures 
usually located above active regions; they ofte n exhibit significant fine structure. 
The moving Type IV bursts (jRobinsonl . 1978) (IV mA) are emitted from sources 
of metre wave continuum moving outwards at velocities of the order of 100- 1000 
kms -1 ; their spectrum is often featureless and sometimes lasts more than ten 
minutes. A number of these are believed to originate wit hin the dense substruc- 
tures (such as erupting prominences) within the CMEs ( Klein and MouradianL 
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2002t iBastian et all , 120011 ). Others appear following Type II bursts and are 
possibl y caused by ene rgetic electrons produced in the wake of the Type II 
shock. iRobinsonl (|1985h calls them flare continua II (FCII) based on the close 
temporal association with the Type II and classifies them in the stationary Type 
IV (IV mB) family with frequency drift. 

Broad-band (Af/f « 1) and narrow-band (Af/f < 0.1) pulsations are, at times, 
considered part of the Type IV family being interpreted as MHD oscillations 
in coronal loops modulating the continuum. Alternative interpretations include 
limit cycles of nonlinear wave-particle interactions in coronal loops or quasi- 
periodic particle acceleration episodes during the magnetic reconnection in a 
large-scale current sheet. On the dynamic spectra, the pulsations exhibit quasi- 
periodic behavior (w 0.1-10 Hz) and durations from seconds to minutes; some 
gradually drift to lower frequency in the course of the event (e.g. reviews by 
Bend . l2003t iNindos and Aurassl . 120071 ). 



Another type of radio burst which has been associated with magnetic recon- 
nection in solar flares is the short (< 0.1 second), narrow-band (some MHz) 
emissions known as "spikes". They appear in broad-band clusters lasting sec- 
onds to about a minute. Metric spikes at 250—500 MHz are quite common 
near the starting frequency of metri c Typ e III bursts and are thus linked with 



the acceleration of electron beams ( Bend . 2003 : Kliem et all . 120031 ). They are 
considered as the signature s of accelerated particles at a highly fragmented, 
primary energy release site ([Nindos and Aurassl . I2007T ) . At the same time they 
pinpoint the location of the reconnection and acceleration site. 

A small number of solar active phenomena has been classified as "extreme 
events"; in these events, characteristics such as flare intensity, CME speed and 
solar-energetic particles (SEP) flux are orders of magnitude larger than the rest 
( Crosby . 20091) . Solar-energetic proton events may result in short duration sharp 
increase in the count rate of ground based cosmic ray detectors and are dubbed 
ground level enhancements (GLEs). 

The extreme event of 20 January 2005 in AR 720 (N14° W61°), included 
a fast CME, an X7.1/2B (6:36-7:26 UT) flare, a white-light flare, hard X- 
rays, 7-rays (up to 300 MeV), and radio bursts in microwave to decametric 
freque ncy. It also initia t ed SEP and GLE, indicating a prompt proton ac celer- 
ation (|Kuznetsov et all 120061 120081 : iGrechnev 6^.120081: Isimnettl . l2007h with 
the hardest energy spectrum of solar cycle 23 ( Bazilevskaval [2009h . The GLE, in 
particular, was characterized as the most intense in half a century and am ong the 
most intense in observation history (|Plainaki efail. 120071 iBelovl . 120051 ). second 
only to the GLE of 23 February 1956 (jVashenvukl . 120051 ). 

This report on the 20 January 2005 extreme event is structured as follows: In 
section [5] we describe the instrumentation and data set used in our study; this 
is supplemented with an algorithm for the detection and statistical analysis of 
quasi-linear structures embedded in complex dynamic spectra in the appendix. 
The observational results are presented in Section[3j they are discussed in Section 
U and compared with models of magnetic reconnection and particle accelera- 
tion within the framework of the CME-induced reconnection. Conclusions are 
presented in Section [3] 



SPsub5B_04arXiv.tex; 5 April 2011; 2:27; p. 3 



4 



C. Bouratzis et al. 



2. Instrumentation and Data Processing 
2.1. Instruments and Data Sets 

The Appareil de Routine pour le Traitement et 1' Enregistrement Magnetique de 
1' Information Spectral (Artemis-IV) solar radio— spectrograph(htt p/ /web.cc.uoa.gr/ artemis/)! 
operating at Thermopylae Greec e (38°49' N, 22°41' E), since 1996 (|Caroubalos et aZj .| 
2001 : Kontogeorgos et all . [20 06') consists of a 7- metre parabolic antenna covering 
the metric range and a dipole antenna covering the decametric range. Two re- 
ceivers operate in parallel, there are a sweep frequency analyzer (ASG), covering 
the 650-20 MHz range in 630 data channels at ten samples sec" 1 and a high 
sensitivity multi-channel acousto-optical analyzer (SAO), which covers the 270- 
450 MHz range in 128 channels at 100 samples sec -1 . The broad-band, medium 
time resolution recordings of the ASG are used for the detection and analysis of 
radio emission from the base of the corona to two Rq , while the narrow-band, 
high time resolution SAO recordings are mostly used in the analysis of the fine 
temporal and spectral structures. 

In this wo rk we used d a ta fro m the Hiraiso Radiospectrograph (Hiras) (25- 
2500 MHz. 



Kondo et all I1995T ) to complement the Artemis-IV range in the 



high-frequency domain. 

The Wind/ Waves experiment ( Bougeret et al , 1995) includes the RAD2 re- 
ceiver which covers the critical frequency range from 13.825 MHz (Rq) to 1.075 
MHz (20 Rq). It is complemented by the RAD1 receiver (frequency range 
1040-20 kHz), which allows us to track the evolution of the radio sources to 
1 AU. 

The dynamic spectra were supplemented with recordings in th e 2-35 GHz 
range of the Nobeyama Radio Polarimetre (NoRP) ( Nakaiima et al , 19851 ). The 
event was also observed on HXR and 7-rays in the 0.03 -300 MeV r ange from 



the Solar Neutron and G amma (SONG)/CORO NAS-F (iKuznetsovL I2004D : we 
used RHESS I data from (|Grechnev efali l2008l ) . The SXR data were obtained 
from GOES (http/ /w ww. sel.noaa.go v/ftpmenu /indices) and the kinematics of the 
CME came from iGopalswamv et all ( 2005 ). 

The times of e nergetic-par t icle release from the low corona were d etermined 
from the reports of lSaid (|2005l ) ; iGrechnev et all (|2008h ; ISimnetd (|2007t) ; iKurt et aL| 
(|2010al fbh. 

Thus, with Wind/ Waves observing the Interplanetary Space and HiRAS with 
Artemis-IV the high-frequency counterpart of the event (corona and lower corona)| 
we obtain a 2000 MHz— 20 kHz frequency coverage. This, combined with the 
SXR, HXR and 7-ray recordings, makes a multi-frequency and multi-instrument 
study of all different aspects of the active phenomena possible, the injection of 
energetic particles included, comprising the 20 January 2005 event. 



2.2. Coronal Density-Height Model Selection 

As plasma emission depends on electron density, which in turn may be converted 
to coronal height using density models, we may calculate estimates for the radio 
source heights and speeds from the dynamic spectra. The establishment of a 
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Figure 1.: Combined HiRAS, Artemis-IV (ASG), and Wind/ Waves dynamic spectrum of the event of 20 January 2005, 06:15- 
07:30 UT and 2000-2 MHz; the CME trajectory using the Newkirk model (see text) for the height to frequency conversion 
(green); the GOES SXR flux (blue); the SONG 40-100 KeV channel (red); the two Type II/FCII combinations (II(1)/FCII(1) 
and II(2)/FCII(2)) ar e in cyan f rames ; det ails are pre sented in Figure |U The electron and proton release times as reported by 
Grechnev et all (|2008h : ISimnettl (l2007h and[Saii (|2005h (Tabled]) are annotated with arrows under the plot. 
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ARTEMIS-IV ASG 2005-20-01 
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Figure 2. Artemis- IV/ASG dynamic spectrum (linear frequency scale 300-40 MHz) of fun- 
damental-harmonic bands of the two metric (11(1) and 11(2)) and the decametric (11(4)) 
shocks 





Figure 3. Narrow-band Type III and spikes high resolution (30 ms) dynamic spectrum (a) 
Artcmis-IV/SAO and the corresponding NoRPTeX microwave enhancement (frequencies 35.0 
(green), 9.4 (blue), 3.75 (red) and 2.2 (white) GHz). This radio-activity in the 06:36-06:44 UT 
interval marks the onset of the 20 January 2005 event, (b) Artemis IV differential spectrum 
(SAO) in the 06:42-06:44 UT interval with 10 msec resolution. On the left a group of spikes 
at 06:42:20 UT, marked by the frame, on the right details of the pulsating patch, (c) Evolution 
of average (logarithmic) frequency drift rate (df/fdt)of the marked spike cluster in the period 
06:42:20-06:42:40 UT; peaks appear at 0.10, ± 0.06, -0.11 sec" 1 . jTsitsipis et am2006l . 120071 , 
also Appendix). 
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Table 1. Overview of the Event 20 January 2005 and associated activity. 



Event 


UT 


Characteristics 


Remarks 


CMh Litt On 
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Gopalswamv et al. (2005) 


Start 


06:36 


3. 1 o-oo.U LrllZ 


INOriJr 


SXK Start 


06:36 


A T> ^70A \T-| AD T T T (~* 1 O 

AR 720: N14 W61 


GOES 12 


IV Start 


06:36 


2000-100 MHz 


Artemis-IV, Hiras 


Narrow-band 


06:39 


Turnover ks200 MHz 


Artemis- IV 


J and U 




Drifting to RT15 MHz 


SAO 


Ha Start 


06:41 


Two ribbon. AR 720 


LEAR 


HXR increase 


06:42 


ai 100 keV 


fMasson et al. 20091 


Spikes Start 


06:42:40 


450-270 MHz 


SAO 


Continuum Patch 


06:43 


460-350(MHz) 
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(F-H) 


to 06:43:45 


910-690(MHz) 


and Hiras 


FCII(l) Start 


06:43 


650 MHz, 


Artemis-IV 






340 Km sec -1 




7-ray increase 


06:43-06:44 


« 780 KeV- 10 MeV 


(Kurt et al, 2010b) 


Microwave Peak 


06:45 


35 GHz 


NoRP 


7-ray increase 


06:45-06:46 


« 60-100 MeV 


SONG/CORONAS-F 








(Kurt et al, 2010b) 


11(4) Start 


06:45:16 


41 (F)-82 (H) MHz 


Artemis-IV 


HXR peak 


06:45:30 


RilOO keV 


RHESSI 


III(l) GG 


06:45:39 




Artemis-IV, Hiras 








J TT/'' J / TT7 

and Wind/ Waves 


11(1) Start 


I \ f ' A p" OA 

06:45:39 


7b{b )-150(H) MHz 


Artemis-IV, Hiras 






120 Km sec 
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fid- A ft Dfi- A 1 
UD.4D-UD.4 ( 


PSlZf — ZZo Rev 


vj-rcennev et at. (zuuoj 
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~1.7 Gcv 
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2B 
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11(1) End 
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UD.Oz 


i on i\/fu„ 


Artemis-IV 


New set of Flare 


06:53-06:57 




TRACE 


Kernels activated 






(cf. Grechnev et al, 2008) 


CME First 


06:54 


PA 293° 


LASCO 


appearance on C2 




3242 Km sec -1 , 


Gopalswamv_et_o/ i (2005) 






or 2075 Km/sec 


Grechnev et al (2008) 


FCII(2) Start 


06:55 


650 MHz, 


Artemis-IV 






380 Km sec" 1 




HXR peak 


06:55 




RHESSI 


(secondary) 






(Grechnev et al, 2008) 


Electron release 


06:55 




Kurt et al. f2010b) 


Microwave peak 


06:55 


9.4 GHz 


NoRP 
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1L(Z ) otart 


UD.OD 
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11(2) End 


06:58 
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11(4) End 


07:02 


36 (F)-72 (H) MHz 


Artemis-IV 


SXR Peak 


07:01-07:12 
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GOES 12 


FCII(2) End 


07:03 


120 MHz 


Artemis-IV 
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07:15 


14 MHz 


Wind/ Waves 






750-4690 Km sec" 1 


Pohiolainen et al. (2007) 


SXR End 


07:26 




GOES 12 


Ha End 


08:54 




LEAR 
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Figure 4.: Artemis— IV (ASG) dynamic spectra of the Type II/FCII bursts (see cyan frames on figure [TJ (a) II(1)/FCII(1) 
06:40-06:50 UT: The electron and proton release times (06:46-06:47 UT, Table [T]) are marked with arrows under the plot. We have 
also included SONG/CORONAS-F normalized flux at 40-100 KeV (blue), 0.775-2.0 MeV (yellow) and 60-100 Mev (red), (b) 
Details of the Type 11(1) shock in the 110-300 MHz range and the 06:45-06:50 UT time interval; it is enclosed by the box in (a), 
(c) II(2)/FCII(2) 06:52-07:02 UT: The electron release time (06:55 UT, Table [1]) is marked with an arrow. We have also included 
SONG/CORONAS-F normalized flux at 40-100 Kev (blue), (d) Details of the Type 11(2) shock in the 110-250 MHz range and 
the 06:55-06:59 UT time interval; it is enclosed by the box in (c). 
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Figure 5. Artemis IV Spectra of II(1)/FCII(1), bidirectional Type III and reverse-drift Type 
Ill-like bursts; (a) ASG dynamic spectrum (06:43-06:47 UT) (b) ASG differential spectrum, (c) 
evolution of average (logarithmic) frequency drift rate (df/fd) of the lowerpart of the differential 
spectrum (06:44-06:46 UT and 450-300 MHz) (Appendix |Aj with peaks at -0.10 sec^ 1 (20 
Mm sec -1 outbound) and 0.83 sec -1 (« 120 Mm sec -1 reverse drift); we have assumed radial 
propagation and a twofold Newkirk density-height coronal model (d) SAO differential spectrum 
06:44:30-06:45:05 UT. 
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frequency of observation-coronal height and frequency drift rate-radial speed 
correspondence is affected by ambiguities introduced by the variation of the 
ambient medium properties. These may be the result of burst exciter propagation 
within the undi sturbed plasma, oyer-dense or under-dense structure or CME 



after-flows (e.g. iPohiolainen et all l2007t iPohiolainen. Hori. and Sakurail 12008 



for a detailed discussion on mo del selection). 



The twofold Newkirk model (|Newkirkl . 1 1961 ) is adequate over active regions 
at heights greater than 1 .05 R p, (frequencies lower than 300 MHz) (figure 4 in 
Wild. Smerd. and Wei ss. 1963), and it is used throughout this report. 



3. Association of Active Phenomena with Radio Bursts 
3.1. Overview 

The complex radio bursts of the e xtreme event of 20 January 2005 have, attracted 
considerable attention (see e.g. |Pohiolainen et all , 120071 : iMasson ~et oD, 120091: 



Grechnev et all l2008t iMavromichalakiL 12009. ). 

The metric-radio-burst complex includes two distinct metric-decametric Type 
II (6:39-06:58 UT), in close succession, each associated in time with a flare- 
continuum- like burst (FCII), and Type III activity both isolated and in large 
groups (III GG). These are superposed on an extended (06:36 UT to well past 
08:00 UT) Type IV continuum with rich fine structure (mostly fibers and pul- 
sations). In the decametric frequency range, a stationary shock (with almost no 
drift) was appears in the 40 (F) and 80 (H) MHz bands from 06:45:12 UT to 
07:02 UT. 

The radio bursts are accompanied by GOES/SXR and HXR, 7-ray emission 
in the 0.03-300 MeV range from SONG/CORONAS-F, and R HESSI; the high- 
energy channels, exhibit a sharp rise of intensity a t 06:45:30 UT ([Kuznetsov et aLl .| 
20081 (Grechnev et all , 120081 IMasson et all . l2009t ). 



Small time-scale narrow-band features were recorded at high resolution by 
Artemis-IV/SAO, as part of the Type IV and FCII fine structures. The latter 
provide evidence of magnetic reconnection and electron acceleration in AR 720 
(N14° W61°). 

An overview of the time sequence of the observations used in this report is 
presented in Table [TJ In Figure [T] we present a combined dynamic spectrum of 
the radio bursts; we have included the CME trajectory, the GOES SXR and 
SONG 40-100 KeV light-curve, and the times of energetic-particle injections for 
comparison. Details of the Type II events (11(1), H(2) and 11(4)) are presented 
in Figures [2] and |H 

In the following subsections, we have partitioned our description of the event, 
and the corresponding radio spectra, for convenience. We present, in subsection 
3.21 the low-frequency part of the dynamic spectrum as it applies to the hecto- 
metric Type II event (11(3)). The description of the radio bursts during the early 
stages of the flare is in subsection 13.31 In the next subsections. 13. 3.11 and l3~3~2l 
we describe the complex II-III GG-FCII bursts and the decametric stationary 
Type II in Section 13.41 Lastly we present, briefly, the Type IV continuum in 
Subsection 13.51 
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3.2. CME-flare Onset and DH Type II Shock 

The event of 20 January 2005 starts with the CME takeoff at 06:33 UT which 
accelerates to high sky-p lane speed at about 06: 54 UT (3242 Km sec" 1 after 
Gopalswamv et all (j2005T) or 2075 Km sec" 1 from iGrechnev et all (|2008t )). The 



GOES 12 SXR flare onset at 06:36 UT reaches X7.1 intensity at 07:01-07:12 
UT (Figure [I| ; its Ha counterpart is a 2B flare from Learmonth (LEAR) in 
06:41-08:54 with a peak at 06:46 UT. At, 06:36 UT the NoRP 3.75-35 GHz flux 
increase starts; this indicates that high-energy electrons were present from the 
very beginning of the event. 

A hectometric Type II event ( 11(3) in Figure [J) appear s at 07:15 UT at 14 



MHz on the Wind/ Waves spectra: IPohiolainen et al. ( 2007 ) calculated the Type 



11(3) speed from the Wind/ Waves spectra to be in the 750-4690 Km sec 1 range 
(depending on the selection of the coronal-density model). 

3.3. The Start of Decimetric-Metric Radio Activity 

The radio activity commences at the same time with the onset of the SXR flare 
and the 3.75-35 GHz enhancement, at 06:36 UT (Figure [3]), with a broad-band 
Type IV continuum which covers, almost, the combined frequency range (2000- 
20 MHz) of Artemis-IV and Hiras (Figured]); its duration extends well past 08:00 
UT (Section 133}. 

A little later (06:39 UT), superposed on the continuum, groups of type J and U 
bursts at decimetric-metric wavelengths appear. Their turnover frequencies drift 
slowly towards the lower part of the spectrum, forming a moving front (Figure 
|3^,). The type J and U bursts are accompanied by groups of narrow-ba nd spikes. 
After the onset of the 100 kcV HXR at 06:40 UT (|Masson et a/.l . l2009l) the spike 



rate increases (06:42:20-06:42:40 UT, Figure^). The spikes exhibit a frequency 
drift rate in the range -100 to 100 MHz sec -1 with distinct peaks at ±0.06 sec -1 , 
±0.10 sec -1 corresponding to exciter speeds in the 10-16 Mm sec -1 range (Figure 
[5k). The Type III(J) bursts and the spikes, mentioned above, possibly mark 
an additional acceleration episode corresponding to a small peak on the NoRP 
recordings. At 06:43:00-06:43:45 UT a continuum patch appears on the dynamic 
spectra in the frequency range 1000-760 MHz (H) and 475-350 MHz (F) (Figure 
m . The continuum patch exhibits a pulsating structure (Figure EU) with periodic 
behaviour (0.32 Hz: main FFT peak, 10 Hz: secondary FFT peak) and coincides 
with the onset of the flare continuum FCII(l). 

3.3.1. The first Type II-FCII-III GG Complex Burst 

A close examination of FCII(l) (Figure [5]) reveals extended groups of Type III 
bursts and reverse-drift Type III bursts sharing a common onset frequency which 
corresponds to the separation frequency of the bidirectional Type III bursts. The 
onset frequency drifts towards lower frequencies. The high-frequency part of the 
FCII(l) (270-450 MHz) was recorded by the Artemis-IV/SAO with high time res- 
ol ution (100 sec -1 ) . Thi s part was examined using the methodology introduced 
by iTsitsipis et all (|2006l 12007 ) which provides statistics on the frequency drift 
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rate based on a 2D FFT of the dynamic spectrum ( Appendix lAf. The logarithmic 
frequency drift rates peak at -0.10 sec -1 and 0.17 sec -1 MHz sec -1 (speeds 15-20 
Mm sec -1 ). A third and more pronounced peak at 0.83 sec -1 is characteristic 
of the reverse-drift Type Ills. For the twofold Newkirk model the exciter speed 
is estimated to ~ 120 Mm sec -1 . 

The Type II (11(1), in FigureH^, and|3}D) appeared on the Artemis-IV record- 
ings at ?» 06:45:40 UT and 145 MHz ending at 06:49 UT and 125 MHz. The 
radial shock speed, as calculated from the twofold Newkirk coronal model, was 
found quite small (« 120 Km sec -1 ). 

The 11(1), FCII(l), and III(l) combination of bursts appears to be associated 
with the rise phase of the SXR flare (Figured]). The intense Type III GG (III(l) 
in Figure[T]) at 06:45:39 UT spans the whole range of the combined Hiras- Artemis 
IV- Wind/ Waves spectra, and it coincides in time with the particle release and 
the HXR-7-microwave peaks as specified in the following paragraph. 

The rise in microwave and 100 keV HXR flux (HXR enhancement starts at 
06:40 UT, peak at w 06:45:30 UT) overlaps in time with the drift of (FCII(l) 
«340 Km sec -1 ) to lower frequencies (6:44-6:46 UT) as it tends to converge 
with 11(1), on the dynamic spectra, near 200 MHz (k, 0.26 Rq, Figures Q] and 
Hk). This FCII(1)-II(1) convergence coincides with a sharp increase of 7-ray 
flux around 6:45-6:46 UT. The latter originates from pion decay, as protons 
with energies in excess of 300 MeV int eract with nuclei i n the ambient corona 
the p ion decay time being w 06:45 UT ( Kuznetsov et al . 20081 : iGrechnev et al 



2008 j) . The pion de c ay em ission supports the path and time-of-flight calculations 



of IGrechnev et al. (|2008l) which provide an estimate of the proton release in 
the 06:46-06:47 UT interval. The electron release is estimated within the same 
interval as the pion decay, and the HXR bremsstrahlung light curves present 
similar profiles. 

3.3.2. The Second Type II-FCII-III GG Complex Burst 

Another Type II occurence (11(2) in Figures [TJ [4J; and HJi starts at 06:56 and 
200 MHz drifting at w 543 km sec -1 . It is accompanied by the FCII(2) event 
(starting 06:55 UT) which consists mainly of a group of reverse-drift Type III 
events with a total drift towards low frequencies corresponding to a radial exciter 
speed ~ 380 Km sec" 1 . It, a lso, coincides with a secondary maximum of the 
HXR ( Grechnev et al . , 2008. their Figure 5) associated with microwave flux 



enhancement, after 6:53 UT; to these correspond an electron release at «06:55 
UT and the second DH Type III GG (111(2)) starting, on the dynamic spectrum, 
from the Type 11(2) lane. This Type II event originated, probably, from a flare 
blast wave. 

The complex II(2)/FCII(2) burst is associated with the maximum of the 
SXR flux but also with a weak HXR, compared to the previous bursts, and 
a microwave peak. 

3.4. The Stationary Decametric Type II 

The Stationary Decametric Type II event starts at 06:45:16 UT at 41 MHz 
(Fundamental), exhibiting a fundamental-harmonic structure; it is labeled 11(4) 
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on Figures [T] and [21 The harmonic structure is not very well defined as it overlaps, 
somewhat, with the fundamental of 11(2) and it coincides, also, in frequency with 
the Artemis-IV data gap. It ends at about 07:02 UT after the secondary energy 
release and the injection of the second Type III group (111(2)). A number of 
Type III bursts appear to originate from the Type II band; they merge into the 
two III GG groups (III(l) and 111(2)). 

3.5. The Type IV Continuum 

The Type IV burst becomes more pronounced after II(2)/FCII(2); it exhibits 
rich fine structure too which includes broad-band pulsations, intermediate drift 
bursts (fibers) and a variety of narrow-band bursts. 

Between FCII(2) and the stationary Type IV a cluster of narrow-band struc- 
tures, mostly Type I II, was recorded by Ar temis-IV/SAO at 300-400 MHz. A 
little after 07:00 UT. iGrechnev et all (2008) observed an expanding arcade, on 



GOES/SXI and EUV (SOHO/EIT) images, which may be the source of the IV 
event. 



4. Summary and Discussion 

We have presented a radio signature analysis of the solar eruptive event of 20 
January 2005 and we have studied the association of the radio bursts with the 
energy release and the active phenomena comprising this event. 

4.1. The Dynamic Spectrum of 20 January 2005 

The two different Type II bursts (11(1) and 11(2)) exhibited significantly different 
drift rates and were well separated in time and each was associated with an HXR 
peak. This suggests distinct shocks. 

The first of them, 11(1), was preceded by metric-decimetric drifting structures 
such as isolated Type III(J) and III(U) bursts, spikes and narrow-band Type 
III events (mostly reverse-slope); these can be identified as the signature of 
reconnection and electron acceleration ep isodes above expanding soft X-ray loops 
(jKlassen. Pohiolainen. and Kleinl . l2003h . Furthermore, the J-U bursts, at the 



beginning of the metric radio event, trace the initial closed magnetic structure. 

Both, 11(1) and 11(2), are accompanied by a flare continuum(FCII(l) and 
FCII(2) respectively) which, on close examination, consists of bidirectional and 
reverse-drift Type III elements; these continua also trace electrons accelerated 
in reconnection sites. The drift rate of the FCII suggests outward movement of 
the reconnection. 

The large Type III groups, III(l) and IH(2), appear to start, in part, from 
the Type II bands (11(1) and 11(2) respectively) and in part from the decametric 
shock 11(4) as discussed in Section 13.41 A few start at higher frequencies; this 
indicates shock acceleration of a significant part of the electrons, exciting the 
Type III events. Prior to the III GGs the Type III(J) or III(U) bursts trace 
energetic electrons which are confined and do not access interplanetary space. 
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Figure 6. Comparison of the standard CME-Flare model with the combined Hiras, 
Artemis-IV (ASG) a nd Wind/ Wa ves dynamic spectrum. On the left we present the CME-in- 
duced reconnection (Forbes, 2003) supplemented with the reconnec tion outflow jets and the 
corresponding termination shocks ( Aurass, Vrsnak, and Mann, 2002) where the upward shock 
appears after the CME. An additional Type II event originating from the flare-loop expansion 
moves sideways. On the right we present the dynamic spectrum resulting from this process. 



The continuation of III(l) and 111(2) into hectometric frequencies implies an 
opening of initially closed magnetic lines and the escap e of en ergetic particles 
into interplanetary space. This agrees with Masson et al. ( 2009t ). as they report, 
before III(l) , a number of acceleration episodes in HXR of which the first few do 
not have DH radio signatures. Within this context, the III(l) and HI(2) onsets 
are consistent with the energetic-particle release times. 



4.2. Acceleration and Release of Energetic Particles within the Framework of 
the Standard Flare-CME model 



The t hree main ca t egorie s of particle acceleration processes (e.g. lAnastasiadisl 



2002 ; Aschwanden , 20081 for a review) in solar flares and CMEs (electric DC 



field acceleration at magnetic X-points, stochastic acceleration by turbulence, 
and MHD shock acceleration in reconnection outflows) are all present in the 
extreme event of 20 January 2005. The DC-field acceleration appears in re- 
connecting current sheet driven by the rising CME. The other two processes 
operate in the reconnection outflow termi nation shock accelerating the highly 
energetic, electron and proton populations ( Mann. Aurass. and Warmuthl . 120061 : 
Mann. Warmuth. and Aurassl 2009t Warmuth. Mann, and Aurassl 20091) . 



SPsub5B_04arXiv.tex; 5 April 2011; 2:27; p. 14 



Radio Observations of the 20 January 2005 X-class flare 



15 



A close association, in time, of the electron injection at «06:45:40 UT, the 
HXR primary peak, and the start of the III GG burst (III(l)) was established 
in Section I3~3TT1 The convergence of II(1)-FCII(1) spectra (06:45-06:47 UT), on 
the other hand, coincid es with the estimated time of proton escape at «06:45:30 
( Grechnev et q/J . 12008) as stated in Section |4~T1 This temporal association seems 
to indicate that 11(1) may be a termination shock formed as the reconnection 
front, FCII(l), rises; on th e other hand the expand ing flare loops may piston 
drive 11(1) as a flare blast ( Vrsnak and Cliver , 2008). In both cases this shock 
formation may be accompanied by particle acceleration and release. What weighs 
more in favour of 11(1) as a flare-associated Type II event is the appearance, in 
the decametric frequency range, of 11(4) whose spectrum appears to be more 
typical of a zero-drift termination shock; this also accelerates a significant part 
of the Type III bursts in the groups 111(1 ) and IH(2). 

In an attempt to introduce a model interpreting the observational data and 
the evolu tion of the event , we fir st ly invoke the C ME driven reconnection as pro - 
posed by iForbes and Linl (|2000l) . iForbesl (|2003l) . and IPriest and Forbesl (|2002h : 
the rising CME stretches and deforms the coronal magnetic-field lines resulting in 
a vertical current sheet, where magnetic reconnection commences. Since X-type 
reconnection is directly associa ted with counterstr e aming electron beams appear- 



ing a s bidirectional Type Ills (jAschwanden et all 11995c lAschwanden and Benz 



19971 ) the FCII(l) is interpreted as a moving reconnection front within the 
fr amework of this mechanism (Figur e [5]) . 

iKliem. Karlickv. and Benzl {2000) interpret pulsating structures such as the 
continuum patch marking the start of FCII(l) in terms of quasi-periodic particle 
acceleration episodes. These originate from a magnetic reconnection in large- 
scale current sheets and are characterized by repeated formation and coalescence 
of magnetic islands, feeding a continuously growing plasmoid which drifts out- 
wards in the corona. In this case, however, there are no further observational 
signatures of a rising plasmoid. 

The slowly drifti ng burst 11(4), may be a reconnec tion outflow term i nation 
shock as proposed bv lAurass. Vrsnak. and Mann ( 2002 ) and Aurass et al. (2003).| 
The only difference is that they envisage a shock formed below the X-reconnection| 
as the outflow encounters the post-flare loop tops; our data, however, suggest 
that the shock is above the reconnection as the Type 11(4) is in the decametric 
range and always is at lower frequencies than the FCII(l). Therefore we con- 
clude that the shock may be formed from the upper outflow (Figure [6]). The 
reconnection outflow termination shock was found to be capable of accelerating 
the, highly energetic, electron and proton popu l ations indicated in the SONG 
recordings (Mann. Aurass. and Warmuthl. 2006 : Mann. Warmuth. and Aurassl 
2009j; Warmuth. Mann, and Aurassl 20091 ). The presence of 11(1) complicates 
the spectral form somewhat; however it may originate at the flanks of the rising 
flare loops, hence it has a rather small speed. 

The difficulties of the detection of the reconnectio n outflow termination shock 
on the top of the post-flare loops are attributed by lAurass. Vrsnak. and Mann 
(|2002l ) to the confusion on the dynamic spectra due to other radio bursts, 
especially during the impulsive phase of the flare. 
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On the le ft of Figure [5] we propose a model including both CME-induced 
reco nnection (Forbes . 2003}) and term ination shocks from reconnection outflow 
jets (jAurass. Vrsnak. and Mannl . 120021) : the combined mechanisms may interpret 
the dynamic spectrum which we have included on the right of the figure for 
comparison. At this point we note that the II(1)/FCII(1)/III(1) complex burst 
was preceded by Type III( J) and Type III(U) bursts, spikes, and the start of the 
Type IV continuum. This indicates that energetic electrons were present before 
the onset of the main part of the energy release in the flare, probably originating 
in the early stages of the reconnection process. 

During the second release of electr ons, manifest as 1 11(2) and associated with 
the secondary HXR peak reported bv lGrechnev et all (|2008l ). it is expected that 
the acceleration region has moved higher as this HXR peak is quite small. The 
Type II burst (11(2)) drifts at «540 Km sec -1 and, being far from the CME, may 
be interpreted as a flare blast at this point. The whole Type Il-flare continuum 
complex (II(2)-FCII( 2)) follows the igniti on of the, so far inactive, flare ker- 
nels (Section 13X2)) . iGrechnev et all (|2008h studied the TRACE 1600 A images 
during the time interval 6:52:41-6:57:30 UT, marked four pairs of kernels and 
computed intensity time profiles of their brightness. As the brightness for two 
pairs of kernels decreases it increases on the other two; this transition appears 
at about 6:53-6:54 UT. We note, at this point, that the decametric shock, 11(4), 
is still present until 07:02 UT, which is almost coincident with the end of the 
II(2)/FCII(2) burst. This relatively l ong duration o f the termination shock radio 
burst has already been reported, by lAurassl ( 20031 ). 



5. Conclusions 

The present study was based on high cadence data from the Artemis-IV solar 
radio-spectrograph in the 650-20 MHz range, extended, as regards the frequency 
range, with Hiras (25-2500 MHz) and Wind/ Waves (13.825-1.075 MHz, 1040- 
20 kHz) recordings and supplemented with NoRP (2-35 GHz), SONG (0.03- 
100 MeV), and GOES data; they were compared with the energetic-particle 
release times reported in the literature. This made possible a multi-frequency 
and multi-instrument study of different aspects of the active phenomena, in the 
20 January 2005 event. Our analysis was focused on the Radio Perspective of the 
event, as radio observations provide a rich collection of diverse diagnostics in the 
study of the shock- formation processes, the particle acceleration, and the mass 
ejection from the Sun. The results were interpreted in terms of the standard 
CME-flare model which applies to eruptive flares and invokes the CME-induced 
reconnection scenario. 

A number of particle accelerators were identified, and their relative impor- 
tance in the manifestation of the event was tested against their radio signatures. 
Firstly the CME-associated shock acceleration seems to be of little importance. 
The reconnection in the wake of the CME and the reconnection outflow termi- 
nation shocks (both above and below the reconnection) appear to be the main 
contributors of energetic electrons and protons. 
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The results are consistent with the standard CME-flare model with the re- 
connection outflow termination shocks and the reconnecting current sheet in the 
wake of the CME as sources of energetic particles. 



Appendix 

A. Detection and Statistical Analysis of Quasi-Linear Structures 
embedded in Complex Dynamic Spectra 

In the study of solar radio bursts, various types of structures often appear on 
the dynamic spectra intertwined, overlapping, disturbed by ter restrial signals , 



and e mbed ded in a continuum background. A methodology (jTsitsipis et al. 



20061 l2007h for the detection and statistical analysis of linear structures has 



been developed and is presented here in brief. The continuous background is 
eliminated by means of high pass filtering; this suppresses slowly varying com- 
ponents enhancing fine structure with fast temporal variation. The terrestrial 
interference, appears in the form of lines parallel to the temporal axis and may 
be suppressed with directional filtering. 

In the resulting clean dynamic spectrum, all frequency channels are normal- 
ized to the same minimum and maximum; this eliminates variations of brightness 
along the frequency axis which may interfere with the subsequent processing 
algorithms. 

For the statistical analysis of quasi-linear structures representing Type III 
events, fibers, pulsations etc., on the dynamic spectra we introduce the concept 
of energy density as a function of angle which measures the signal power of a 
2D image (a dynamic spectrum in this case) along a certain direction. In terms 
of the two-dimensional Fourier transform of the dynamic spectrum: 

oo oo 

F(£i,h)= J I f(x,y)-e-i 2 ^ +y ^dxdy (Al) 

— oo — OO 

The energy density as a function of angle (angular energy density) is defined as: 

OO 

S(6) = 2j\F(£,e)\ 2 Zdt (A2) 
o 

This exhibits local maxima at the points where 9 equals the slopes of linear or 
quasi-linear segments within the image, yet it is not affected by their position 
within it. These peaks of energy per direction on the image indicate, in a statisti- 
cal sense, the dominant slopes of linear or quasi-linear segments within the image, 
or, in this case, the peaks of frequency drift rates which may be interpreted as 
the dominant exciter velocities once a coronal density-height model has been 
adopted. An example of the usage of this method appears in subsection 13.31 and 
in Figure [5] (panels (c) and (e)), where we present the statistics of the frequency 
drift rates. 
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